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ABSTRACT 

With the discovery of Kuiper Belt binaries that have wide separations and 
roughly equal masses new theories were proposed to explain their formation. Two 
formation scenarios were suggested by Goldreich and collaborators: In the first, 
dynamical friction that is generated by a sea of small bodies enables a transient 
binary to become bound (L 2 s mechanism); in the second, a transient binary gets 
bound by an encounter with a third body (L 3 mechanism). 

We show that these different binary formation scenarios leave their own unique 
signatures in the relative abundance of prograde to retrograde binary orbits. This 
signature is due to stable retrograde orbits that exist much further out in the Hill 
sphere than prograde orbits. It provides an excellent opportunity to distinguish 
between the different binary formation scenarios observationally. 

We predict that if binary formation proceeded while sub-Hill velocities pre- 
vailed, the vast majority of all comparable mass ratio binaries have retrograde 
orbits. This dominance of retrograde binary orbits is a result of binary forma- 
tion via the L 2 s mechanism, or any other mechanism that dissipates energy in a 
smooth and gradual manner. For super-Hill velocities binary formation proceeds 
via the L 3 mechanism which produces a roughly equal number of prograde and 
retrograde binaries. These predictions assume that subsequent orbital evolution 
due to dynamical friction and dynamical stirring of the Kuiper belt did not alter 
the sense of the binary orbit after formation. 

Subject headings: Kuiper Belt — planets and satellites: formation 
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INTRODUCTION 



The detection of comparable mass binaries with wide separations in the Kuiper Belt 



called for new theories explaining their formation (e.g. 


Weidenschilline 


2002; 


Goldreich et al. 


2002; 


Funato et al. 


2004; 


Astakhov et al. 


2005: 


Lee et al. 


2007). Their existence cannot be 



explained with a formation scenario involving a c ollision and tidal evolution, as has been pro- 
posed for the formation of the Moon and Charon ( Hartmann fc Davisll975 ; Cameron fc Ward 
19761 ; lMcKinnonlll989l ). since it cannot provi de the current angular momentum of the binary 
system. In a formation scenario proposed by IWeidenschiilingj (120021 ) two Kuiper Belt objects 
(KBOs) collide with each other inside the Hill sphere of a third. However, in the Kuiper Belt, 
gravitational scattering between the two intruders is about 100 timea^ more common than 
a collision. The r efore, three body gravitational deflection (L 3 mechanism), as proposed by 
Goldreich et al.l (120021 ). should dominate the binary formation o yer such collisional scen ario. 
A second binary formation scenario that has been suggested by I Goldreich et al.l (120021 ) con- 
sists of the formation of a transient binary that gets bound with the aid of dynamical friction 
from a sea of small bodies. We call this the L 2 s mechanism. In the formation scenario of 



Astakhov et al.l (120051 ) and iLee et al.l (120071 ) the existence of long lived transient binaries 
that spend a long time in their mutual Hill sphere, ne ar a periodic orbit, i s responsible for 
the creation of Kuiper Belt binaries (KBBs). Finally, iFunato et al.l (120041 ) proposed a bi- 
nary formation mechanism that involves a collision between two large KBOs. This collision 
creates a small moon that is replaced in an exchange reaction by a massive body with high 
eccentricity and large semi-major axis. 

In this paper, we show that the L 2 s and L 3 mechanism leave unique signatures in the rel- 
ative abundance of prograde to retrograde binary orbits. The L 2 s m echanism dominates over 
the L 3 mechanism for sub-Hill velocities (jSchlichting &; Saril 120081 ). We argue that binaries 
that form from dynamically cold KBOs by the L 2 s mechanism have retrograde orbits. This is 
due to the existence of stable retrograde binary orbits with modified Jacobi constants similar 
to that of unbound KBOs on circular orbits that have impact parameters that correspond 
to distances o 
bits exist (e.g. 



' closes 


t approach of less than the Hill radius. No equivalent prograde or 


Henon 


1970; 


Innanen 


1979; 


Zhang & Innanen 


1988; 


Hamilton & Burns 


1991 



Hamilton fc Krivovl 119971 ). Since dynamical friction only gradually increases the modified 
Jacobi constant (for a binary this corresponds to gradually increasing the absolute value of 
the binding energy), all binaries that form via the L 2 s mechanism, or any other mechanism 
that dissipates energy in a smooth and gradual manner, will start with modified Jacobi 



1 For this estimate we used a ~ 10 4 and assumed that the velocity dispersion of the KBOs at the time 
of binary formation is less than their Hill velocity, see §2 for details 
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constants close to that of unbound KBOs that penetrate the Hill sphere and hence have 
retrograde orbits. For super- Hill KBO velocities, on ly the L 3 mechanism can form tight 
binaries that tend to survive (jSchlichting fc Sari 2008 ). The fact that retrograde orbits are 



stable for larger semi-major axes is no longer of importance since only tight binaries are 
saved from break up. This, therefore, leads to the formation of a roughly equal number of 
prograde and retrograde binaries for super-Hill KBO velocities. 

Our paper is structured as follows: In §2 we outline our assumptions, explain our choice 
of parameters and define variables that will be used throughout this paper. We calculate the 
ratio of prograde to retrograde binary orbits for the L 2 s and L 3 mechanism and predict the 
relative abundance of prograde to retrograde orbits for sub-Hill and super-Hill KBO velocities 
in §3. We compare our predictions with observations in §4. Discussion and conclusions follow 
in 85. 



DEFINITIONS AND ASSUMPTIONS 



The Hill radius denotes the distance from a body at which the tidal forces due to the Sun 
and the gravitational force due to the body, both acting on a test particle, are in equilibrium. 
It is given by 



vti\ + m 2 
3M 



1/3 



where m\ and m 2 are the masses of the two KBOs, a is their semi-major axis around the 
Sun and the mass of the Sun. Our definition of the Hill radius differs from that used by 
Schlichting fc Saril (120081 ) since we include the combined mass of both KBOs here. We chose 
to do so since it will make comparisons with works by other authors easier. 



We use the 'two-group approximation' (IGoldreich et al.l 120021 . 120041 ) which consists of 
the identification of two groups of objects, small ones, that contain most of the total mass 
with surface mass density a, and large ones, that contain only a small fraction of the total 
mass with surface mass density S C a. We assume a ~ 0.3g cm -2 which is the extrapola- 



tion of the minimum-mass solar neb ular ([ Hayashi 



Estimates from Kuiper Belt surveys (ITruiil. 



o et al. 



1981) to a heliocentric distance of 40AU. 



200 ll ; iTruiillo & Brownll2003l : IPetit et al. 



20081 : iFraser et alJ 12008 : iFuentes fc Holmanll2008r ) yield S ~ 3 x 10~ 4 g cm~ 2 for KBOs with 
radii of R ~ 100 km. We use this value of S, assuming that £ during the formation of 
KBBs was the same as it is now. Our choice for £ and a is also consistent with results from 



numerical coagulation simulations by lKenyon fc Luul (119991 ) . 



Large bodies grow by the accretion of small bodies. Large KBOs viscously stir the small 
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bodies, increasing the small bodies' velocity dispersion u. As a result u grows on the same 
timescale as R provided that mutual collisions among the small bodies are not yet important. 
In this case, u is given by 



where a = R/R H ~ 10~ 4 at 40AU (IGoldreich et all 12002 ). v H is the Hill velocity of the 



large bodies which is given by Vh = QRh where Q is the orbital frequency around the sun. 
The velocity v of large KBOs increases due to mutual viscous stirring, but is damped by 
dynamical friction from a sea of small bodies such that v < u. Balancing the stirring and 
damping rates of v and substituting for u from equation ([2D, we find 

— ~a~ 2 f-^) -0.1. (3) 
v H V / 

For our choice of parameters, we have sub-Hill KBO velocities during the epoch of formation 
of bodies with R ~ 100km. We therefore focus our work on the shear-dominated velocity 
regime (v <C vh)- However, we discuss how our results would be modified if v 3> vh- 



PROGRADE VERSUS RETROGRADE BINARY ORBITS 



3.1. Sub-Hill velocities: w<C% 



The disk of KBOs is effectively two-dimensional in t he shear-dominated yelocity regime 
(v -C vh), since the growth of inclinations is suppressed (IWetherill fc Stewartlll993l ; iRafikov 
20031 ; IGoldreich et al.l I2004J ) . We therefore restrict our calculations for the shear-dominated 
velocity regime to two dimensions. Since we are intereste d in close encounters among the 
KBOs, their interaction is w ell described by Hill's equations (jHilllll878l : IGoldreich fc Tremaine 
19801 ; iHenon fc Petitl Il986l ). In Hill coordinates the equations of motion of the two KBOs 
can be decomposed into their center of mass motion and their relative motion with respect 
to one another. The modified Jacobi constant is exactly conserved in the Hill formalism, but 
the Hill formalism itself is an approximation to the general three body problem. It assumes 
that the masses of body 1 and 2 (in our case the two KBOs) are much less t han that of the 



Sun. We us e th e stan dard Hill coordinate system and reference frame as in IHenon fc Petit 
(119861 ) and lldal (119901 ). In this rotating frame the direction of the x-axis is given by the 



line connecting the Sun and the center of mass of the two KBOs such that the positive x 
direction is pointing away from the Sun. The y-axis is perpendicular to the x-axis pointing 
in the direction of motion of the KBOs' center of mass around the Sun. In Hill coordinates 
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the modified Jacobi constant is 



J c = 3x 2 + 



6 



( X 2 + y 2y/2 



■ 2 -2 

-a; -y 



(4) 



where x and y correspon d to the relative sep aration between the two KBOs in the x and 
y direction respectively (IHenon fc Petitl Il986l ). Length has been scaled by Rh and time 
by Q^ 1 . In Hill coordinates the Lagrangian points L\ and L 2 are located at (—1,0) and 
(+1, 0) respectively, where we define L\ as the Lagrangian point located between the KBO 
and the Sun. Their modified Jacobi constants are Jc{L\) = Jc(L 2 ) = 9. From equation 
H] we can see that tight binaries with small separations have Jc S> 9. We call a binary 
orbit prograde if its angular momentum about the binary center of mass, as viewed in the 
non rotating frame, is in the same direction as the orbital angular momentum of the binary 
around the Sun. If the binary angular momentum is in the opposite direction to the orbital 
angular momentum of the binary around the Sun, the orbit is called retrograde. Several 
authors recognized t hat planar retrograde orbit s are stable for larger semi-major axes than 
prograde orbits (e.g. lHenon|ll970l ; Ilnnanenlll979l ; IZhang fc Innanenlll988l ; Hamilton fc Burns 



199ll ; lHamilton fc Krivovlll997l ). A prograde binary with an initially circular orbit becomes 
unbound for a > 0A9Rh where a is the initial semi-major axis of the mutual binary orbit 



(lHamilton &: Burnslll99ll ). This implies that prograde orbits with modified Jacobi constants 
less than that of the Lagrangian points L\ and L 2 are unbound. In contrast to the prograde 
case, there exist stable retrograde binary orbits with J c < Jc{L\) = Jc(L 2 ) = 9. This result 
is also shown in Figures [H Figure [I] shows histograms of Jc for prograde and retrograde 
binaries that formed by L 3 mechanism from KBOs with initially circular orbits around the 
Sun. In the reminder of this paper we discuss the stability of prograde and retrograde orbits 
in terms of Jc and not semi-major axis since the latter is not well defined (i.e. it is not 
a constant of motion) for wide orbits with a ~ Rh- The modified Jacobi constant for two 
KBOs that approach each other from infinity is 



J 



c 



3x z 



x 



y 



3 -* 

4 



(5) 



where b is the initial separation between the two KBOs in the x direction and e is the relative 
eccentricity in Hill units given by |ei — where ei and e2 are the eccentricity vectors of 
body 1 and body 2 respectively. Only KBOs with b ranging from 1.7Rh to 2.5Rh penetrate 
each others Hill sphere if started on circular orbits. From equation ([5]) we have therefore that 
only KBOs with 2.2 < Jc < 4.7 have a distance of closest approach of Rh or less provided 
that they started on circular orbits around the Sun. 
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3.1.1. L 2 s Mechanism 

In the L 2 s mechanism KBBs form from transient binaries that become bound with the 
aid of dynamical friction from a sea of small bodies. This dynamical friction provides a 
gentle force that damps the random velocity of large KBOs. For typical parameters, the 
dynamical friction force only extracts a small fraction of energy over an orbital timescale. 
Therefore, KBBs that form via the L 2 s mechanism, or any other mechanism that dissipates 
energy gradually, have initially modified Jacobi constants similar to that of the unbound 
KBOs that penetrate within the Hill sphere. As mentioned above, for KBOs that started 
on circular orbits around the Sun this corresponds to 2.2 < Jc < 4.7. However, only stable 
retrograde orbits exist with Jc < 9. This implies that all KBB that form this way must 
have retrograde orbits since no stable prograde orbits exist for Jc < 9. Once a binary is 
formed, dynamical friction increases the modified Jacobi constant and the absolute value of 
the binary binding energy. We confirm that all binaries that form from KBOs on initially 
circular orbits around the Sun via the L 2 S mechanism are retrograde by numerical integrations 
that are presented below. 

Since it is not feasible to examine the interactions with each small body individually, 
their net effect is modeled by an averaged force which acts to damp the large KBOs' non- 
circular velocity around the Sun. We parameterize the strength of the damping by a di- 
mensionless quantity D defined as the fractional decrease in non-circular velocity due to 
dynamical friction over a time 

D-^(-Yor 2 ^^-a 2 (^-Y . (6) 



pR \vh J pR \Vh 

The first expression is simply an estimate of dynamical friction by a sea of small bodies 
assuming u > vh- The second expression describes the mutual excitation among the large 
KBOs for v <C vh- v achieves a quasi-steady state on a timescale shorter than at which R 
grows since only a subset of the deflected bodies are accreted. The stirring among the large 
KBOs can th erefore be equat e d to the damping due to dynamical friction (for a detailed 



derivation see iGoldreich et al.l (120041 ) ). 



Since the growth of inclinations is suppresse d in the shear-dominated v elocity regime 



the disk of KBOs is effectively two-dimensional (IWetherill fc Stewartl Il993l ; iRafikovl 12003 



Goldreich et al.l 120041 ) . We therefore restrict this calculation to two dimensions. In Hill 
coordinates the relative motion of two equal mass KBOs, including the dynamical friction 
term, is governed by 

x - 2y - 3x = y 2~\ — 2W2 _ D± ( 7 ) 



(x 2 + y 2 fl 2 

Length has been scaled by Rh and time by VL~ l . Equations (EJ) and (jHJ) are integrated for 
different values of D and impact parameters ranging from 1.7Rh to 2.5Rh with equal step 
size. Impact parameters outside this range result in a distance of closest approach between 
the two KBOs of more than Rh- 

For D = 0.01, we performed 20000 integrations. About two percent of these integrations 
resulted in the formation of a binary. Figure [2] shows three examples of the evolution of the 
specific angular momentum and Jc of binary formation events from our integrations for 
D = 0.01. In addition, we performed integrations for values of D ranging from 0.1 to 0.0004 
and find that, just like in the D = 0.01 case, only retrograde binaries form. We define h as 
the specific angular momentum of the binary in the non rotating frame. It can be written 
as h = xy — yx + x 2 + y 2 and is related to the total binary orbital angular momentum, 
L, by h = (1/mi + l/m 2 )L. The time t — corresponds to the time at which y — 
if the relative KBO velocity is solely due to the Keplerian sheer (i.e. ignoring the actual 
gravitational interaction between the bodies). The evolution of h and Jq is shown until 
the binary separation has decreased to O.IRh or less. Binaries with separations of O.IRh 
or less are sufficiently tight that perturbations from the Sun are too weak to flip the sign 
of the angular momentum. As expected from our discussion above, the angular momenta 
of the binaries are negative corresponding to retrograde binary orbits. In fact all binaries 
that form via the L 2 s mechanism in our numerical integrations display retrograde orbits. 
Dynamical friction shrinks the binary separation. As a result, the magnitude of the binary 
angular momenta decreases with time. The right hand side of Figure [2] shows the evolution 
of the modified Jacobi constant. Newly formed binaries initially have a modified Jacobi 
constant < 9 which is possible only for retrograde binaries. Dynamical friction shrinks the 
semi-major axes of the binaries which leads to an increase of Jc with time while keeping the 
sense of rotation, i.e. the sign of h, fixed. Eventually the modified Jacobi constant grows 
to values above Jc(Ll) = Jc(L2) = 9. For J c > 9 prograde orbits can exist; however all 
binaries that formed with the aid of dynamical friction started out with Jc < 9 for which 
only retrograde orbits are stable. Therefore, all KBB that form via the L 2 s mechanism, 
or any other mechanism that gradually removes energy from transient binaries, orbit each 
other in the retrograde sense since otherwise they would not be able to form in the first place. 
Figure [3] shows the evolution of h and Jc as a function of time for KBO encounters that 
did not lead to the formation of a binary. These examples show that KBOs encounter each 
other and leave each other with positive angular momenta. This is a result of the Keplerian 
sheer and follows from the definition of h. 

We have assumed here that all KBOs are initially on circular orbits around the Sun 
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and have shown that this leads to the formation of exclusively retrograde binaries in the L 2 s 
mechanism. If, however, the velocity dispersion of the KBOs is sufficiently large, such that e 
is of the order of the Hill eccentricity, bigger impact parameters allow the KBOs to penetrate 
each others Hill sphere. In this case, there now exist KBOs that have an initial Jq just a 
little below 9 (see equation [5]) in which case only a small change in Jc is sufficient for the 
formation of retrograde and prograde binaries. Therfore, prograde binaries can form with 
the aid of dynamical friction provided that the velocity dispersion of the KBOs is about Vh- 

Our prediction for the sense of the binary orbit relies on the assumption that dynamical 
friction does not alter the sense of the binary orbit in the subsequent binary evolution. 
Although we have shown in our simulations that for our dynamical friction model this is 
indeed the case, it might be that the actual behavior of dynamical friction differs from the 
model implemented here. 



3.1.2. L 3 Mechanism 



A transient binary forms when two large KBOs penetrate each other's Hill sphere. This 
transient binary must lose energy in order to become gravitationally bound. In the L 3 
mechanism the excess energy is carried away by an encounter with a third massive body. 
This encounter can provide a significant change in energy which corresponds to a considerable 
change in Jc- The modified Jacobi constants of KBBs that form via the L 3 mechanism are 
therefore not constraint to values similar to that of their initial Jc', their orbits can therefore 
be both prograde and retrograde. We show that this is indeed the case with numerical 
integrations discussed below and determine the ratio of prograde to retrograde orbits for 
binary formation via the L 3 mechanism. 

Our calculation is performed in the shear-dominated velocity regime in two dimensions. 
As initial condition, we assume that all bodies are on circular orbits. We modify Hill's 



equations (IHillll 1 8781 : 1 Goldreich fc Tremaindll980l ; lPetit fc Henonlll986l ) to include three equal 
mass bodies besides the Sun. The equations of motion, with length scaled by Rh and time 
by for body 1 are given by 



Xi - 2y x - 3xi 



3(xi - x 2 ) 



3(xi - x 3 ) 



2((*i - x 2 y + ( Vl - y 2 ) 2 )3/ 2 2((xx - x 3 y + ( Vl - y 3 ) 2 ) 3/2 



3(2/1 - 



3(yi - ys) 



2((xi - x 2 f + ( Vl - y 2 ) 2 ) 3 / 2 2(( Xl - x 3 ) 2 + (yi - y 3 ) 2 ) 3/2 ' 



(9) 



(10) 



The subscripts 1, 2 and 3 label the x- and y-coordinates of KBO 1, 2, and 3 respectively. 
Similar equations of motion can be obtained for bodies 2 and 3. Resulting binary orbits 
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are calculated by n umer ically integrating the equations of motion. We refer the reader to 



Schlichting fc Saril (120081 ) for the exact details of these calculations. 



Figure [T] shows histograms of the modified Jacobi constant of prograde and retrograde 
binaries that formed via the L 3 mechanism. Both histograms are normalized to unity. As 
discussed above, we indeed find that prograde orbits only exist for Jc > 9. The stability 
of retrograde orbits extends below Jc = 9 down to Jc ~ — 10. It therefore includes the 
range 2.2 < Jc < 4.7 with orbits that penetrate the Hill sphere from circular heliocentric 
orbits. Unlike the L 2 s mechanism, the L 3 mechanism does produce retrograde and prograde 
binaries for w<C%. We find that 65% of all binary orbits are retrograde and 35% prograde 
(see Figure H]). Here, we only considered binary formation from three equal mass bodies that 
started on initially circular orbits around the Sun. We therefore caution, that the ratio of 
prograde to retrograde orbits due to the L 3 mechanism might differ for other mass ratios 
and velocity dispersions. 



3.1.3. The Ratio of Retrograde to Prograde Orbits 



Schlichting fc Saril (120081 ) have shown that for sub-Hill KBO velocities the ratio of the 



L to L s binary formation rate is 



FR Li v , . 

^_i- = 0.05— . 11 



Therefore, for sub-Hill KBO velocities, binaries in the Kuiper Belt form primarily due to 
dynamical friction. For our estimate of (v/vh) ~ 0.1, we have that FR l a / FR L 2 S ~ 0.005, 
in which case ~ 0.5% of all binaries form directly by the L 3 mechanism. Since prograde 
binaries can only form via the L 3 mechanism, they make up a negligible fraction of the total 
binaries. Below we discuss how a somewhat larger fraction of prograde binaries can arise 
due to exchange reactions with unbound KBOs. 

Once a binary is formed its semi-major axis shirks due to dynamical friction provided 
by a sea of small bodies. Dynamical friction decreases the orbit of a KBB that has an orbital 
velocity % at a rate 

n shrink ~ Dtt ~ ^-a- 2 [ — J (12) 
pR \v H J 

where we assume that vb < u. Exchange reactions or binary break up by passing KBOs 
occurs at a rate given by 

^■exchange ~ ^ a_2 ^ (jfj ' 
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The ratio of these two rates is given by 

T^shrink ( 



"^exchange ^ V 



(14) 



where v <C v h and Vb > Vh- Break up or exchange reactions are most likely for wide binaries, 
in which case vb ~ vh since vb increases as the semi- major axis of the mutual binary orbit 
decreases. Therefore we have from equation (fill) that H s hrink /^-exchange ~ (vh/ v ) ~ 10 for 
our estimate of (v/vh) ~ 0.1. This implies that ~ 10% of all binaries that formed will suffer 
an exchange reaction or break up. We performed numerical integrations of binary break up 
and exchange reactions to obtain a more accurate estimate and find that only about 3% of the 
binaries suffer an exchange reaction and/or break up. Our order of magnitude calculation, 
therefore, slightly over estimates the number of binaries that experience an exchange reaction 
and/or break up. Moreover, only a fraction of the these binaries will end up as binaries with 
prograde orbits. In conclusion, we predict that the vast majority (> 97%) of comparable 
mass ratio binaries will have retrograde orbits if KBO velocities of v < O.Ivh prevailed 
during binary formation. This prediction assumes that subsequent orbital evolution due to 
dynamical friction did not alter the sense of the binary orbit after formation. 



3.2. Super-Hill Velocity: v > v H 



There is some uncertainty in what the actual values of a and £ were during binary 
formation. For a few times larger value of £ with a unchanged, we enter the regime in which 
v exceeds the Hill velocity (this can be seen from equation d3J)). We discuss here briefly how 
this would affect the ratio of prograde to retrograde binary orbits. 



Schlichting fc Saril (120081 ) have shown that, for v ^> Vh, only binaries that form with a 
binary separation of R cr u = Rh(vh/v) 2 or less tend to be saved from break up. The L 2 s 
mechanism fails in creating binaries with separations ~ Rcru or less since dynamical friction 
is not able to dissipate sufficient energy for tight binaries to form. Therefore, the L 2 s mech- 
anism is not important if KBOs have super-Hill velocities. Tight binaries (with separations 
less than < R cr it), can form via the L 3 mechanism. However in this case, the binary forma- 
tion cross section is s ignificantly reduced with respect to the sub-Hill velocity regime (see 



Noll et al.l (120071 ) and ISchlichting fc Saril (120081 ) for details). The fact that retrograde orbits 



are stable for larger semi-major axes is no longer of importance since only tight binaries tend 
to survive. We therefore predict that a roughly equal number of prograde and retrograde 
binaries form if super Hill velocities prevail. This prediction is supported by Figure HI Figure 
H] shows the ratio of retrograde binaries with a modified Jacobi constant of J™ m or larger 
to the total number of binaries that formed via the L 3 mechanism for v <C vh- When all 
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binaries are included we find that about 2/3 have retrograde orbits. More retrograde than 
prograde binaries form because retrograde binary orbits are stable further out in the Hill 
sphere than prograde ones. As J™ m increases the fraction of retrograde binaries decreases 
reaching a minimum of about 1/3 for J™ m ~ 9. This may be due to the Keplerian sheer 
which increases the duration of a prograde encounter between unbound KBOs compared to 
a retrograde encounter. The fraction of prograde and retrograde binaries becomes compa- 
rable for J™ n 3> 9 because for such binaries neither the Keplerian sheer nor the increased 
stability of retrograde orbits are important. This is the relevant regime for binaries that 
form for v ^> v h since these large modified Jacobi constants correspond to tight binaries, 
which are the only binaries that are saved from break up if super-Hill velocities prevail. 



4. COMPARISON WITH OBSERVATIONS 



To date more than a dozen KBBs have well determined orbits (e.g. iNoll et al.l 120071 ). 
Unfortunately due to projection effects, the prograde and retrograde orbital solutions of 
the KBBs are nearly degenerate. This degeneracy can usually only be broken after several 
years once the viewing angle of the KBBs has changed sufficiently Very recently, after the 
submission of our original manuscript, two groups reported unique orbital so lut ions for KBBs 



Typho n-Echidna flGrundv et all 120081 ) and 2001 Q W™ flPetit et al. 



20081 ). iGrundv et al. 



(120081 ) find a prograde orbit for Typhon-Echidna and IPetit et al.l (120081 ) report a retrograde 
orbit for 2001 QW^- 2001 QW322 has such a large binary separation that, even in the current 
Kuiper belt, it experiences significant dynamical interactions with other large KBOs. It is 
early to draw conclusions for the whole binary population, but if a comparable number of 
retrograde and prograde binaries is found it would imply that KBBs formed from a dispersion 
dominated KBO disk, which would also be consistent with observed binary inclinations. 
Dispersion dominated KBO velocities would imply that the value of H/a was larger during 
binary formation than what we used in equation [31 However, the velocity dispersion during 
binary formation cannot have exceeded Vm significantly sin c e the binary formation ti mescales 
would otherwise become excessively long (INoll et al.l 120071 ; ISchlichting fc Sarill2008l ). 



5. DISCUSSION AND CONCLUSIONS 

The relative abundance of prograde to retrograde orbits enables us to differentiate be- 
tween various proposed binary formation scenarios observationally. 

We predict that the vast majority (> 97 %) of comparable mass ratio binaries will 
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have retrograde orbits if KBO velocities of < O.Ivh prevailed during their formation. This 
dominance of retrograde over prograde binary orbits is due to the fact that for sub-Hill 
velocities binaries form primarily via the L 2 s mechanism instead of L 3 mechanism. Since 
dynamical friction only gradually increases the modified Jacobi constant, all binaries that 
form via the L 2 s mechanism, or any other mechanism that dissipates little energy over an 
orbital timescale, will start with modified Jacobi constants close to that of unbound KBOs. 
Only stable retrograde orbits exist for modified Jacobi constants similar to that of KBOs 
with initially circular orbits around the Sun that penetrate inside the Hill sphere. Therefore, 
KBBs have retrograde orbits provided that they form from dynamically cold KBOs via the 
L 2 s mechanism. 

As the KBO velocities approach vh the preference of retrograde orbits decreases. Fur- 
ther, we predict a comparable number of prograde and retrograde binaries form for super-Hill 
KBO velocities. This is because only th e L 3 mechanism can for m tight binaries that tend 



to survive if super Hill velocities prevail (jSchlichting fc Sarill2008l ). The fact that retrograde 



orbits are stable for larger semi-major axes is no longer of importance since only tight bi- 
naries tend to survive. This therefore leads to the formation of a roughly equal number of 
prograde and retrograde binaries for super-Hill KBO velocities. 

The analysis pr esented here has also im plications for some of the other proposed binary 



formation scenarios. IWeidenschillingi (120021 ) suggested that KBBs form by a collision among 
two KBOs inside the Hill sphere of a third. Although the L 3 mechanism dominates over such 
a collisional binary formation scenario we briefly discuss our predictions for this collisional 
binary formation mechanism. For sub-Hill velocities more retrograde than prograde binaries 
form because retrograde binary orbits are stable further out in the Hill sphere than prograde 
ones (i.e. the phase space for forming retrograde binaries is larger than that for prograde 
binaries). For super-Hill velocities a comparable number of prograde and retrograde binaries 
form because the fact that retrograde orbits are stable for larger semi- major axes is no longer 
of importance since on ly tight binaries are saved from break up. In the formation scenario of 



Astakhov et al.l (120051 ) the existence of long lived transient binaries that spend a long time 



in their mutual Hill sphere, near a periodic orbit, is responsible for the creation of KBBs. 



Lee et al.l (120071 ) find an excess of prograde over retrograde binaries and suggest that this is 
a signature of their binary formation process. Our work indicates that an excess of prograde 
over retrograde binaries might simply be the result of the velocity regime (i.e.v ~ vh) in 
which the binaries form (see Figure H|) . 

All of the above predictions rely on the assumption that subsequent orbital evolution 
due to dynamical friction and dynamical stirring of the Kuiper belt did not alter the sense 
of the binary orbit. The Kuiper Belt has undergone a phase of dynamical excitation which 
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probably modified the orbital properties of KBBs. A detailed study on how dynamical 
stirring of the Kuiper Belt and dynamical friction affects binary inclinations would be very 
worthwhile to determine whether they were able to reverse the binary orbit from prograde 
to retrograde rotation. 

Some of the numerical calculations presented here were performed on Caltech's Division 
of Geological and Planetary Sciences Dell cluster. R. S. is an Alfred P. Sloan Fellow and a 
Packard Fellow. 
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Fig. 1. — Histogram of modified Jacobi constants, Jc, of prograde and retrograde KBBs that 
formed via three body gravitational deflection, L 3 mechanism, for u < %. Each histogram is 
normalized to unity, but overall retrograde orbits are twice as abundant as prograde orbits. 
Note, prograde binaries exist only for J c > 9 whereas retrograde binaries exist also for 
Jc<9- 
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Fig. 2. — Three examples of KBO encounters in the L 2 s mechanism for v <C v h and D = 0.01 
that result in the formation of a binary. The plots on the left and right hand side show the 
evolution of the specific angular momentum, h, and modified Jacobi constant, Jc, as a 
function of time respectively. The time t — corresponds to the time at which y = 
if the relative KBO velocity is solely due to the Keplerian sheer (i.e. ignoring the actual 
gravitational interaction between the bodies). The evolution of h and Jc is shown until 
the binary separation has decreased to O.IRh or less. These examples show that the sense 
of rotation is practically preserved, h displays large variations right after capture caused 
by solar tides. The most extreme case of angular momentum sign change for bodies that 
form binaries found in our simulations is displayed in the second of the three examples. The 
angular momenta of the binaries are all negative corresponding to retrograde binary orbits. 
In fact all binaries that form via the L 2 s mechanism in our numerical integrations display 
retrograde orbits. Dynamical friction shrinks the binary separation leading to a decrease in 
the magnitude of h and an increase of Jc with time. The modified Jacobi constant of the 
newly formed binaries is smaller than Jc(Ll) = 9 which explains why all their orbits are 
retrograde (see §3.1.1. for details). 
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Fig. 3. — Same as in Figure[3]but for two examples of KBO encounters in the L 2 s mechanism 
for t)<% and D = 0.01 that do not result in the formation of a binary. As a result of the 
Keplerian sheer, KBOs encounter and leave each other with positive h. 
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Fig. 4. — The ratio of retrograde binaries, N ret , with a modified Jacobi constant of J™ m or 
larger to the total number of binaries, N tota i, that formed via the L 3 mechanism for t; < 
For small J™ n , i.e. when all binaries are included, about 2/3 have retrograde orbits. More 
retrograde than prograde binaries form because retrograde binary orbits are stable further 
out in the Hill sphere than prograde ones. As J™ n increases the fraction of retrograde 
binaries decreases reaching a minimum of about 1/3 for J™ n ~ 9. This may be due to 
the Keplerian sheer which increases the duration of a prograde encounter between unbound 
KBOs compared to a retrograde encounter. The fraction of prograde and retrograde binaries 
becomes comparable for J™ n 3> 9 because for such binaries neither the Keplerian sheer nor 
the increased stability of retrograde orbits are important. 



